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ATMOSPHERIC  MEASUREMENTS  NEAR  118  GHz 
WITH  PASSIVE  MICROWAVE  TECHNIQUES 


David  H.  Staelin 
Philip  W.  Rosenkranz 
Alan  Cassel 
David  McDonough 
Paul  Steffes 


1 . 0 Introduction 

Meteorological  measurements  today  are  made  primarily  by 
in  situ  sensors  located  on  the  ground,  ships,  aircraft,  bal- 
loons, and  rockets.  Such  measurements  are  limited  by  eco- 
nomics to  observations  well  separated  in  time  and  space, 
especially  over  oceans  and  in  the  Southern  Hemisphere.  This 
rather  serious  lack  of  global  coverage  has  been  remedied  in 
part  by  satellites  which  photograph  and  measure  the  electro- 
magnetic spectrum  of  the  atmosphere  so  as  to  yield  informa- 
tion about  clouds  and  the  temperature  and  humidity  distribu- 
tions of  the  atmosphere.  All  such  optical  observations  have 
a serious  problem  in  their  inability  to  penetrate  deeply 
most  types  of  clouds.  Although  high  spatial  resolution  has 
enabled  optical  devices  to  sense  between  clouds,  this  tech- 
nique doesn't  work  well  when  the  clouds  are  too  extensive. 
Such  extensive  cloud  cover  is  frequently  found,  for  example, 
near  severe  storms  where  observations  are  of  great  interest. 

Fortunately  most  such  clouds  can  be  penetrated  by 
microwave  radiation,  and  satellite-borne  passive  microwave 
sensors  can  accomplish  many  of  the  objectives  sought  with 
infrared  atmospheric  sounders. 

This  superior  ability  of  microwave  sensors  to  penetrate 
clouds,  particularly  cirrus,  has  motivated  the  development 
of  microwave  spectrometers  for  emote  sensing  of  the  atmo- 
sphere from  space.  On  11  December  1972  the  first  microwave 
atmospheric  temperature  sounder  was  launched  on  the  Nimbus-5 
experimental  earth-observatory  satellite  and  it  has  operated 
continuously  since  then.  This  initial  experiment  has 
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essentially  confirmed  theoretical  predictions  that  a three- 

channel  spectrometer  operating  near  0 . 5-cm  wavelength  can 

yield  atmospheric  temperature  profiles  0-20  km  with  1-3  °K 

rms  accuracy.  The  instrument  has  also  yielded  measurements 

of  atmospheric  liquid  water  and  water  vapor  over  ocean  and 

other  geophysical  parameters.  The  preliminary  results  of 

this  experiment  were  discussed  by  Staelin  et  al . ^ Soviet 

microwave  spectrometers  operating  at  longer  wavelengths 

yielded  information  about  atmospheric  water  and  surface 

2 

parameters  for  several  days  from  the  Cosmos-243  and 
Cosmos-384^  satellites,  which  were  launched  in  1968  and 
1970,  respectively. 

The  initial  satellite  experiments  are  just  beginning  to 
exploit  the  microwave  spectrum  for  passive  microwave  remote 
sensing.  They  have  utilized  only  the  water  vapor  resonance 
at  22.235  GHz  and  the  low  frequency  wings  of  the  60-GHz 
oxygen  absorption  band.  Plans  are  now  being  made  to  use 
such  sensors  on  future  series  of  operational  meteorological 
satellites,  such  as  Tiros-N  and  the  DMSP  Block  5D  satel- 
lites . 

Stronger  resonances  of  water  vapor,  oxygen,  ozone  and 
other  molecules  occur  at  shorter  wavelengths,  1-5  mm,  and 
more  precise  measurements  of  clouds  are  also  feasible. 
Shorter  wavelengths  furthermore  permit  higher  angular  reso- 
lution with  smaller  antennas,  which  can  be  very  important 
for  geosynchronous  satellites. 

The  118  GHz  resonance  of  oxygen  is  of  interest  for 
five  reasons:  1)  the  diameter  of  a 118  GHz  antenna  need  be 
only  half  that  required  for  a 60  GHz  temperature  sounder  in 
order  to  achieve  some  desired  spatial  resolution, 

2)  simpler  radiometers  can  be  used  at  the  118  GHz  resonance 
of  (>2  than  at  the  60  GHz  complex  because  only  a single 
local  oscillator  is  required,  3)  the  118  GHz  resonance  is 
more  intense  than  those  near  60  GHz  and  has  less  Zeeman 
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splitting;  therefore  it  can  be  used  to  sound  higher  in  the 
atmosphere,  4)  the  different  dependence  on  temperature  and 
pressure  may  offer  small  remote  sensing  advantages,  and 
5)  cloud  profiles  can  be  more  accurately  extracted  from  a 
combination  of  60  and  118  GHz  spectral  data  than  from 
60  GHz  data  alone. 

In  this  report  are  discussed  1)  ground-based  observa- 
tions at  60  and  118  GHz,  2)  the  aircraft  program  at 
118  GHz,  and  3)  the  ability  of  118-GHz  systems  to  contrib- 
ute to  linear  retrievals  of  atmospheric  temperature, 
humidity,  and  cloud  profiles. 

2 . 0 Ground-Based  Observations  near  60  and  118  GHz 

2 . 1 Introduction 

The  ground-based  observation  program  was  designed  to 
address  three  questions:  1)  Is  the  118  GHz  line  shape 
predicted  by  present  theory?  2)  Is  the  60  GHz  line  shape 
predicted  by  present  theory?  and  3)  Do  these  two  bands 
respond  to  clouds  as  expected? 

2 . 2 The  60-GHz  Observing  System 

The  60  GHz  system  (AAFE  Microwave  Temperature  Sounder) 
is  a 12-channel  Dicke-switched  radiometer.  It  is  similar  to 
the  NASA  Nimbus-6  SCAMS  instrument  in  that  it  has  a dual 
antenna  system  which  scans  by  rotating  a hyperboloid  reflec- 
tor about  the  axis  of  the  feed  horn.  Figure  2.1  is  a block 
diagram  illustrating  the  double-sideband  superheterodyne 
scheme.  Note  that  an  orthomode  transducer  is  associated 
with  each  antenna  feed  horn,  and  that  channels  1-7  involve 
a second  superheterodyne  stage.  Dicke-switching  is  accom- 
plished with  ferrite  switches.  A particularly  stable  local 
oscillator  (necessary  to  straddle  the  two  oxygen  absorption 
lines)  is  provided  by  crystal-controlled  oscillators  with 
phase-locked  loops.  All  other  hardware  is  standard. 
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The  r.f.  components  are  mounted  on  a cart  which  can  be 
moved  to  any  desired  viewing  position.  The  components  on 
the  cart  receive  power  from  and  return  data  to  a large  equip- 
ment rack  which  contains  hardware  to  process  the  digital  data 
and  to  control  the  mode  of  operation.  On  top  of  the  cart  are 
two  scanning  reflectors,  situated  so  as  to  provide  viewing 
angles  as  shown  in  Fig.  2.2.  Note  the  two  calibration  tar- 
gets, which  are  equipped  with  monitoring  thermistors,  and 
are  mounted  on  the  cart.  Actually,  each  reflector  has  a 
separate  ambient  target.  The  single  heated  target  is  a 
large  hollow  absorbing  wedge  controlled  to  ^ 350  °K  tempera- 
ture. The  cart  and  equipment  rack  are  connected  by  five 
large  cables  which  allow  for  moving  the  cart  on  the  observ- 
ing platform  while  keeping  the  remaining  equipment  indoors. 

Table  2.1  gives  some  specifications  for  the  system  as 

4 

reported  by  its  developers  at  the  Jet  Propulsion  Laboratory. 

Independent  measurements  of  AT  were  carried  out  first  by 
5 rms  ~ 

Badian,  and  then  again  in  August,  1977.  The  1977  measure- 
ments are  included  in  Table  2.1.  These  measurements  were 
performed  by  placing  a large  piece  of  Eccosorb  ^ absorbent 
material  in  front  of  the  antenna,  producing  a constant 
temperature  (ambient,  303.7  °K)  target.  The  resulting  rms 
temperature  fluctuations  are  a factor  of  2 to  4 higher  than 
those  of  the  functional  specifications.  No  concerted 
efforts  have  been  made  to  find  and  eliminate  the  noise 
source,  because  the  noise  is  adequately  reduced  by  increas- 
ing the  integration  time.  The  noise  apparently  arises  be- 
cause of  degraded  mixer-amplifiers.  For  a Dicke-switched 
radiometer , 

* 2TsySte»,//5TI 

where  TSyStem  is  the  equivalent  system  noise  temperature  (K), 
B is  the  bandwidth  (Hz) , and  is  the  integration  time 
(seconds) . 
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Table  2.1.  AAFE  System  Specifications 


Chan. 


Bandwidth 

0.65 

MHz 

1.0 

MHz 

1.0 

MHz 

2.0 

MHz 

2.5 

MHz 

3.6 

MHz 

10.0 

MHz 

50.0 

MHz 

100 

MHz 

100 

MHz 

100 

MHz 

100 

MHz 

atrms  (ti  = 1 Sec-) 


3.71  K 
2.99  K 
2.99  K 
2.11  K 
1.89  K 
1.57  K 
0.94  K 
0.42  K 
0. 30  K 
0.30  K 
0.30  K 
0.30  K 


9.85  K 
11.22  K 
7.57  K 
5.97  K 

4.19  K 
3.60  K 

2.20  K 
1.01  K 
0.95  K 
0.82  K 
0.73  K 
0.90  K 


Frequency 


± 50  KHz 


± 1 MHz 
± 5 MHz 


2.3  The  118-GHz  System 

The  118-GHz  system  is  also  a Di.cke-switched  double- 
sideband superheterodyne  system  (see  block  diagram. 

Fig.  2.3).  A mechanical  Dicke-chopping  wheel  is  used  with 
a rotating  reflector  in  front  of  the  antenna  feed  horn.  As 
the  reflector  rotates  about  the  axis  of  the  feed  horn,  it 
scans  the  angles  shown  in  Fig.  2.4  and  views  the  two  cali- 
bration targets,  which  are  monitored  by  thermistors.  At 
present,  the  local  oscillator  is  a Klystron  tube,  but  pro- 
visions have  been  made  for  the  installation  of  a solid-state 
L.O. 

The  front-end  components  up  to  and  including  the 
detectors  and  75-dB  video  amplifiers  are  mounted  on  a cart. 
Support  equipment  is  rack-mounted,  just  as  in  the  60-GHz 
system,  so  the  two  systems  operate  similarly.  The  118-GHz 
system  also  has  an  ASR-33  teletype  connected  to  the  support 
equipment  rack  for  operator  control  of  data  collection  modes. 
Operation  of  the  118-GHz  system  was  analyzed  by  Paul  G. 
Steffes7  in  the  summer  of  1977. 

2 . 4 Observations  at  60  and  118  GHz 

Three  days  of  observations  have  been  chosen  as  being 
representative  of  the  data  collected  over  several  months. 
Table  2.2  lists  the  dates,  times  of  observations,  and 
weather  conditions  for  the  chosen  data.  Discussion  of  the 
systems'  performances  will  accompany  the  presentation  of 
the  data.  On  all  three  occasions,  the  60  GHz  system  was 
operated  in  its  full  13-position  scanning  mode,  although 
only  the  zenith-looking  data  will  be  analyzed,  and  the 
118  GHz  system  was  operated  at  a viewing  angle  10°  from 
zenith. 

2.4.1  Error  Analysis 

The  60  GHz  system  exhibits  very  quiet  performance. 
Inspection  of  the  raw  data  has  shown  that  all  calibration 


Block  diagram  of 


ZENITH 


Figure  2.4.  118  GHz  system  scan  angles. 


Table  2.2.  Time  and  Weather  for  Observations 


Date 

Time 

Weather 

October  31,  1977 

4:20  pm  E .S . T. 

Clear 

Surface  Temp.  = 284  K 
Dew  Point  = 281  K 

| 

November  15,  1977 

12  noon  E.S.T. 

Variable  Cloudiness 
Surface  Temp.  = 281  K ' 
Dew  Point  = 278  K 

I 

December  18,  1977 

11:40  am  E.S.T. 

Cloudy  (thick) , Light 
Snow. 

Surface  Temp.  = 273  K ' 
Dew  Point  = 272  K 

constants  vary  slowly  compared  to  the  sampling  intervals 

used  to  measure  them.  Small  mismatches  between  predictions 

and  measurements  can  be  accounted  for  by  adjusting  the  so- 

called  nonresonant  linewidth  parameter  (Av  typically  0.58 

MHz  mbar-"*")  used  in  the  temperature  prediction  programs. 

This  parameter  may  be  thought  of  as  a free  variable  in  the 

overlapping  line  model  of  the  60  GHz  band.  Changing  Av  to 

^ 0.48  MHz  mbar  ^ reduces  the  predicted  T for  channel  10 

B 

by  'v  10  K and  reduces  Tfi  for  channels  1-9  and  11  approxi- 
mately half  as  much. 

The  118  GHz  system  unfortunately  did  not  exhibit  quiet 
performance  until  the  end  of  the  program.  The  data  to  be 
presented  will  therefore  have  large  error  brackets. 

Some  of  the  errors  were  systematic.  First,  the  radio- 
sonde measurements  are  made  at  Chatham,  Massachusetts,  which 
is  ^ 75  km  from  our  laboratory.  Consequently,  the  assumed 
Cambridge  temperature  profile  may  be  erroneous.  By  observ- 
ing actual  predictions  for  many  days,  it  can  be  estimated 
that  these  errors  may  be  'v  ±2  K.  Second,  there  is  a 5% 
uncertainty  in  the  best  laboratory  measurements  of  the 
oxygen  absorption  coefficient.  This  translates  to  'v  ±1  K 
potential  errors.  Third,  errors  in  the  frequency  of  the 
local  oscillator  of  'v  0.1  GHz  may  introduce  a +3  K error  in 
Tg.  Fourth,  diurnal  variations  in  the  temperature  profile, 
which  should  not  affect  the  high-altitude  channels,  may  be 
as  much  as  5 K on  channels  1 and  2,  and  2 K on  channel  3. 
Finally,  there  is  evidence  of  a systematic  miscalibration 
of  the  system  involving  the  thermistors  monitoring  the 
calibration  targets.  Small  errors  in  this  calibration  are 
magnified  through  extrapolation  by  as  much  as  a factor  of 
four  when  Tg  'v  160  K.  This  error  could  be  'v  +20  K for 
channel  7 and  *v  +4  for  channel  1.  Figure  2.5  shows  the 
total  error  bracket  for  each  channel. 
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2.4.2  Clear  Air  Results 


Data  were  taken  under  clear  skies  simultaneously  at  60 
and  118  CHz  on  October  31,  1977.  Figures  2.6  and  2.7  show 
the  60  and  118  GHz  results,  respectively.  The  general 
shapes  of  the  curves  agree  qualitatively  with  theory.  The 
largest  contribution  to  uncertainty  in  the  60  GHz  data  is 
receiver  noise  (ATDM_  is  listed  in  Table  2.1)  which  is  re- 
duced  by  integrating  for  52  seconds.  Note  the  general 
features  of  the  Tg  spectrum.  The  coldest  values  of  Tg  occur 
on  channels  for  which  the  atmosphere  is  most  transparent, 
particularly  channel  10.  They  will  be  referred  to  as  "high- 
altitude"  channels.  Data  from  the  118  GHz  system  is  treated 
similarly,  with  higher  numbered  channels  also  sensing  higher 
altitudes.  Channel  1 of  the  118  GHz  system  was  overwhelmed 
by  the  UHF  television  interference  in  the  i.f.  passband. 

This  interference  probably  degraded  channels  6 and  7 also. 
The  striking  difference  between  the  performances  of  the  two 
systems  is  indicated  by  the  error  bars  on  the  118  GHz  data. 

2.4.3  Results  for  Variable  Cloudiness 


On  November  15,  1977,  data  were  collected  with  both  the 
60  and  118  GHz  systems  as  thin  patches  of  clouds  passed  over- 
head. Figure  2.8  shows  that  the  liquid  water  content  of 
these  clouds  did  not  have  a substantial  effect  on  the  mea- 
surement of  Tg  at  frequencies  near  60  GHz.  Correspondingly, 
Fig.  2.9  shows  no  measurable  effect  on  the  118  GHz  data, 
although  calibration  errors  may  have  overwhelmed  a poten- 
tially measurable  effect.  Channels  6 and  7 suffered  inter- 
ference from  television  signals  once  again. 

2.4.4  Results  for  Heavy  Overcast 

Observations  made  on  December  18,  1977,  under  heavy 
overcast  and  some  light  snow  show  distinct  effects  which  may 
be  attributed  to  nonresonant  absorption  by  liquid  water 
(clouds) . Data  from  the  60  GHz  system  (see  Fig.  2.10)  show  a 


j 
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Channel 


O Zenith  Prediction 


Figure  2.6.  60  GHz  Data  Oct.  31,  1977 
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9 X increase  in  Tg  for  channel  10  over  the  predicted  clear- 

sky  value,  and  about  half  as  much  for  intermediate  height 

channels  1-9  and  11.  This  is  in  good  agreement  with  the 

effects  expected  for  uplooking  measurements.  Despite  the 

large  error  bars  for  the  118  GHz  spectrum  (see  Fig.  2.11), 

it  also  shows  a similar  effect.  On  channel  6 (channel  7 

suffered  some  television  interference) , the  increase  in  Tg 

is  between  8 K and  36  K,  probably  much  higher  than  the 

effect  at  60  GHz.  The  values  of  T approach  predictions  for 

lower  altitude  channels,  as  they  should.  If  we  assume  the 

transparent  channels  near  53  and  118  GHz  were  altered  9 K 

1 3 

and  25  K,  respectively,  this  corresponds  to  'v  v ' . This 
frequency  dependence  is  less  than  one  expects  for  droplets 
much  smaller  than  a wavelength,  but  is  quite  reasonable  for 
these  short  wavelengths  and  heavy  clouds. 

3 . 0 Observations  near  118  GHz  from  Aircraft 

3 . 1 System  Description 

The  radiometer  was  completed  in  February,  1977,  and 
was  configured  as  illustrated  in  Fig.  2.3.  It  was  packaged 
in  two  racks  and  a window  panel,  all  of  which  were  provided 
by  the  NASA  Ames  Research  Center.  The  window  panel  con- 
tained the  scanning  antenna,  local  oscillator,  mixer,  and 
i.f.  amplifiers.  See  Figure  3.1.  These  assemblies  were 
then  shipped  to  Moffett  Field,  California,  in  February  for 
installation  on  NASA's  Convair  CV-990  aircraft,  together 
with  other  microwave  and  optical  sensors.  The  radiometer 
was  connected  to  the  NASA/Goddard  Data  System,  the  aircraft 
digital  data  acquisition  system  (ADDAS) , and  the  aircraft 
time-code  generator  (Fig.  3.2). 

The  data  were  collected  and  stored  on  tape  in  blocks  of 
1024  8-bit  words.  Each  1.6  second  integration  produced  two 
numbers  per  channel;  one  for  the  received  power  with  the 
Dicke  chopper-wheel  open  and  one  for  it  closed.  One  such 
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Main  control  panel 
for  118-GHz  radiometer 
showing  the  tape 
recorder,  instrument 
control  unit,  spectrum 
analyzer,  antenna 
control,  teletype, 
power  panel,  and 
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and  shielded  i.f. 
amplifiers.  Hose 
supplied  dry  gas 
to  prevent  conden 
sation  inside 
antenna  assembly. 


Figure  3.1.  Convair-990  Aircraft  Antenna  (viewing 

Installation  of  calibration  load) 

118-GHz  Radiometer.  mounted  in  window. 

Note  rain  shields. 
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Figure  3.2b.  Data  Block  Format 
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block  of  data  required  between  80  and  120  seconds  to  take, 
the  variation  depending  on  the  amount  of  teletype  output 
requested  before  the  next  data  block  could  begin.  The  tele- 
type output  was  controlled  by  front  panel  switches,  as  was 
the  antenna  scan  angle.  The  teletype  and  digital  tape  were 
supplemented  by  a one-channel  strip-chart  recorder  and  an 
equipment  log  book. 

There  were  three  principal  problems  with  the  operation 
of  the  system.  First,  there  was  a recurring  problem  with  an 
i.f.  2.0  GHz)  oscillation  which  saturated  our  video- 
amplifiers on  the  1891  MHz  and  total  power  channels.  Second 
was  a software  problem  complicated  by  the  third  problem, 
intermittent  behavior  of  the  magnetic  tape  drive.  The  net 
effect  of  these  problems  was  to  render  about  half  of  the 
data  unusable. 

3 . 2 Atmospheric  Observations 

Sixteen  data  flights  were  made,  these  are  listed  in 
Table  3.1.  The  accuracy  of  the  aircraft  data  was  limited 
by  the  receiver  thermal  noise  and  by  various  systematic 
errors  as  discussed  earlier.  Linear  combination  of  all 
sources  of  error  yields  3a  uncertainties  of  +4.9,  -1.9  K for 
high-altitude  measurements.  Figures  3.3  and  3.4  present 
comparisons  of  theoretical  and  observed  118-GHz  spectra 
obtained  from  the  Convair  990  data  at  altitudes  of  9.4  and 
5.5  km.  These  figures  appear  in  Paul  Steffes'  master's 
thesis7 (Department  of  Electrical  Engineering  and  Computer 
Science,  September  1977). 

The  data  suggest  that  there  is  no  major  discrepancy 
between  theory  and  experiment  for  absorption  near  118  GHz 
in  the  troposphere  and  that  the  core  of  the  atmospheric  line 
at  low  pressures  may  possibly  be  more  opaque  than  expected 
(see  Fig.  3.3).  Confirmation  of  these  results  should  be  an 
objective  of  future  experiments  of  this  type. 


f 
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Table  3.1.  Flight.  Schedule  Summary- 


Flight  Operators*  Mission  Data  Status 


#1 

3/8/77 

PS,  JS 

Ames  to  Steamboat  Springs, 
Colorado,  to  Ames. 

Manual  calibration  only.  : 

ADDAS  not  recorded. 

#2 

3/11/77 

PR 

Ames  to  NOAA  Ocean  Data 

Buoy  EB-16,  to  Ames. 

Down-looking  calibration 
only.  ADDAS  not  recorded. 

#3 

3/16/77 

PR 

Ames  to  ocean  near  Southern 
California,  to  Ames. 

Same  as  02. 

#4 

3/17/77 

PR 

Ames  to  ocean  near  Southern 
California,  to  Ames. 

Same  as  02. 

#5 

3/18/77 

PR 

Ames  to  ocean,  seeking 
Cirrus,  to  Ames. 

Tape-drive  failure  (no 
data  recorded) . 

#6 

3/21/77 

PR,  AC 

Ames  to  McChord  Air  Base, 
Washington . 

Tape-drive  intermittent, 

3 tapes  used. 

07 

3/22/77 

PR,  AC 

McChord  Air  Base  to  EB-19, 
to  McChord . 

Same  as  02. 

08 

3/23/77 

AC,  JB 

McChord  Air  Base  to  EB-16, 
to  Ames. 

Same  as  02. 

09 

3/24/77 

JB 

Ames  to  ocean  near  Southern 
California,  to  Ames. 

Same  as  02. 

#10 

3/25/77 

JB,  DM 

Ames  to  ocean  near  Southern 
California,  to  Ames. 

Same  as  // 2. 

#11 

3/28/77 

JB,  DM 

Ames  to  ocean  near  Southern 
California,  to  Ames. 

Full  calibration. 

#12 

3/30/77 

JB,  DM 

Ames  to  McChord,  McChord  to 
Elmendorf  Air  Base. 

Same  as  #11.  ADDAS 
recorded  2nd  half. 

#13 

4/1/77 

DM 

Elmendorf  to  Thule  Air  Base. 

ADDAS  recorded.  Full 
calibration. 

#14 

4/2/77 

DM 

Thule  local  flight,  going 
north. 

Same  as  #13. 

#15 

4/4/77 

DM 

Thule  local  flight,  going 
south. 

Same  as  #13. 

#16 

4/5/77 

DM 

Thule  to  Malmstrom  Air  Base, 
Malmstrom  to  Ames. 

Same  as  #13. 

* 

* 

PS 

« Paul 

Stef fes 

JB 

- John 

Barrett 

1 

PR  = Philip  Rosenkranz 
AC  « Alan  Cassel 
DM  * David  McDonough 
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Figure  3.4.  Plotting  of  observed  "brightness  temperatures 
:~rp~ rrr-  r~  and  theoretical  temperatures  derived  from  model 
!_□_ : j.3  versus  channel  number.  • > ; •-» 
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4 . 0 Theoretical  Evaluation  of  Atmospheric  Temperature, 

Humidity,  and  Cloud  Retrievals 

4 . 1 Introduction 

A regression  method  was  used  to  develop  a linear,  least 
squares  fit  to  atmospheric  parameters  based  on  theoretical 
brightness  temperature  observations.  The  ability  of  various 
passive  microwave  sounders  to  support  retrievals  of  these 
parameters  was  evaluated  for  different  climatologies. 

4 . 2 Regression  Technique 

The  regression  technique  used  was  as  follows:  If  we 
let  p be  the  unknown  parameter  vector  (a  column  vector)  and 
d be  the  observed  data  vector,  our  estimate  is  of  the  form 

p*  = Dd 

where  p*  is  the  estimated  value  of  p and  D is  a two-dimen- 
sional matrix  selected  so  as  to  minimize  the  mean  square 
difference  between  p and  p* . To  add  generality,  the  data 
vector  d is  augmented  by  the  inclusion  of  a constant  term 
(in  our  case,  1),  which  has  the  effect  of  adding  constant 
terms  from  the  D matrix  to  the  estimate  p* . Otherwise,  we 
would  be  constrained  to  having  p*  = 0 when  d = 0,  where  0 
denotes  the  zero  vector;  this  is  not  generally  correct  or 
desirable. 

In  order  to  determine  the  matrix  D defined  above,  we 
may  differentiate  the  mean  square  error  with  respect  to 
for  all  i and  j,  and  set  the  result  equal  to  zero.  Thus, 

3I77  E f (p*-p)  (p*-p)  = E (dV'-pS  (Dd-p)  I 

ij  1 ij  > 

= ^ — Etd^Dd  - dtDtp  - pfcDd  + pfcp] 


= E[2djDid  - 2d.p.]  = 
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1 


where  E denotes  the  expected  value  of  its  argument,  and  t 
denotes  the  transpose. 

The  symbol  denotes  the  tth  row  of  the  D matrix. 

Thus , 

E[dpt]  = E[ddtDt]  k CdDfc 

where  denotes  the  data  correlation  matrix.  Hence,  we 
obtain 

D = {C'1  E[dnfc] }t 
as  the  solution  for  D. 

This  approach  to  estimation  is  suboptimum  for  several 
reasons:  1)  actual  statistics  are  not  jointly  Gaussian, 

2)  the  statistics  are  non-stationary  in  space  and  time, 

3)  it  normally  neglects  the  additional  information  content 
arising  from  the  horizontal  correlations  in  the  atmosphere, 
and  4)  the  equations  of  radiating  transfer  are  generally 
non-linear.  The  linear  approach  used  here  does,  however, 
establish  an  upper  bound  on  the  expected  errors. 

4 . 3 Explanation  of  Channel  Sets 

A program  exists  to  obtain  the  theoretically  observable 
brightness  temperatures  from  specified  atmospheres.  This 
program  uses  the  frequencies  designated  in  Table  4.1. 
Channels  1-11  comprise  what  is  referred  to  here  as  the 
118  GHz  set.  Although  weighting  functions  for  these  chan- 
nels were  not  explicitly  calculated,  the  weighting  functions 

O 

for  a similar  set  were  calculated  by  Pettyjohn  and  are 

Q 

presented  in  Figure  4.1.  (See  Staelin,  1969  for  a general 
discussion  of  microwave  weighting  functions.)  The  set  used 
here  differs  from  that  set  slightly  in  that  some  of  the 
channels  have  been  shifted  slightly  to  provide  somewhat  more 
emphasis  on  accurate  soundings  in  the  lower  atmosphere. 

Channel  12  is  the  100  GHz  window  channel,  and  channel 
13  is  the  135  GHz  window  channel.  Channels  14-39  comprise 
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Table  4.1.  Frequency  Designations 


1) 

118.725  GHz 

10) 

115.95  GHz 

20) 

50.30 

2) 

118.6975 

11) 

115.15 

21) 

53.74 

3) 

118.65 

12) 

100.00 

22) 

54.96 

4) 

118.55 

13) 

135.00 

23) 

57.888185 

5) 

118.35 

14) 

182.3725 

24) 

57.722185 

6) 

118.15 

15) 

181.435 

25) 

57.672485 

7) 

117.85 

16) 

180.185 

26) 

57.639'185 

8) 

117.50 

17) 

178.935 

27) 

22.20 

9) 

• 116.75 

18) 

177.06 

28) 

31.60 

19) 

174.56 

(KM) 
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the  183  GHz  channel  set,  which  is  intended  to  sense  water 
vapor.  Channels  20-26  comprise  an  augmented  Tiros-N  Micro- 
wave  Sounder  Unit  (MSU)  set,  and  this  group  is  referred  to 
here  as  the  60  GHz  set.  Channels  27  and  28  comprise  the 
20  GHz  channel  set,  which  consists  of  one  22  GHz  water  vapor 
channel  and  one  31  GHz  window  channel. 

4 . 4 Explanation  of  Atmospheric  Data  Sets 

The  radiosonde  data  sets  which  were  used  to  compute  the 
correlation  matrices  C^  and  E[pdt]  consist  of  400  observa- 
tions in  each  of  three  regions;  arctic  (60°-90°  N/S) , mid- 
latitude (30°-60°  N/S),  and  tropical  (30°  S-30°  N) . The 
radiosondes  were  launched  at  various  locations  within  these 
ranges  and  at  various  times  during  the  day  and  year  over  the 
period  1966-1969.  These  1200  radiosonde  data  records  reside 
on  a single  9-track  tape.  The  arctic  data  set  was  not  used 
in  these  experiments.  The  locations  of  the  stations  used  to 
compile  the  data  sets  are  listed  in  Table  4.2.  Stations  are 
included  in  the  arctic,  midlatitude,  or  tropical  data  sets 
based  on  their  latitude,  as  above,  except  that  station 
number  5 was  categorized  as  a middle  latitude  station. 

For  processing,  only  100  of  the  400  atmospheres  were 
selected  at  a time.  They  were  processed  into  files  of  100 
atmospheres,  which  then  contain  atmospheric  temperatures  and 
computed  water  vapor  content  at  the  following  levels: 

1000  mb,  850,  700,  500,  400,  300,  250,  200,  150,  100,  70, 

50,  30,  10,  and  5 mb.  Cloud  models  were  added  to  some  of 
the  data  sets  as  described  elsewhere  in  this  report.  In 
addition  to  the  observed  radiosonde  parameters,  theoretical 
brightness  temperatures  at  all  28  frequencies  were  calcu- 
lated; that  is,  TB  , TB(jown'  an^  atmospheric  opacity  were 
included  for  each  atmosphere,  for  each  of  28  channels.  In 
addition,  surface  reflectivity  and  temperature  were  calcu- 
lated for  sea  retrievals;  for  land  retrievals,  different 
values  were  used. 
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Table  4.2.  Station  Directory 
Number  Name  and  Location 


1 

Thule,  Greenland 

77N 

69W 

2 

Barking  Sands,  Hawaii 

22N 

160W 

3 

White  Sands,  New  Mexico 

3 2N 

106W 

4 

Primrose  Lake,  Canada 

55N 

HOW 

5 

Fort  Greely,  Alaska 

64N 

146W 

6 

Ascension,  South  Atlantic 

8S 

14W 

7 

Point  Mugu , California 

34N 

119W 

9 

Fort  Churchill,  Canada 

59N 

94W 

10 

West  Geirinish,  Scotland 

57N 

7W 

11 

Wallops  Island,  Virginia 

38N 

7 5W 

12 

Antigua,  West  Indies 

17N 

62W 

13 

Cape  Kennedy,  Florida 

2 8N 

81W 

14 

Eniwetok,  Marshall  Is. 

11N 

162E 

15 

Fort  Sherman,  Canal  Zone 

9N 

8 0W 

18 

Kwajalein,  Marshall  Is. 

9N 

168E 

19 

Natal,  Brazil 

6S 

35W 

21 

Thumba,  India 

9N 

77E 

22 

McMurdo,  Antarctica 

78S 

165E 

23 

Heiss  Island,  USSR 

8 IN 

58E 

24 

Volgograd,  USSR 

49N 

45E 

25 

Mar-Chiquita , Argentina 

38S 

57W 

26 

Molodezhnaya , Antarctica 

68S 

46E 

27 

Chamical,  Argentina 

30S 

6 6W 

41 

Research  Ship  (USSR) 

30-60S , i 

4 . 5 Equation  of  Radiative  Transfer 

Table  4.3  is  a list  of  the  routines  used  in  computing 
the  atmospheric  brightness  temperatures.  References  to  the 
literature  are  also  listed.  The  five  parameters  TB  , 
TBdown'  0PAC>  REFL»  and  TS  represent  upward  brightness 
temperature  above  the  atmosphere,  downward  brightness 
temperature  at  the  surface,  atmospheric  opacity,  surface 
reflectivity,  and  surface  temperature,  respectively.  They 
were  combined  to  give  brightness  temperatures  as  follows: 

TB  = TB  + OPAC*  (REFL*TB.  + (1  - REFL)  *TS). 
up  down 

In  the  case  of  ocean  estimates,  TS  and  REFL  were  computed 
using  the  programs  described  above.  In  the  case  of  land 
estimates,  TS  was  assumed  to  be  the  same  as  the  1000  mb 
temperature  and  REFL  was  computed  as  a Gaussian  random 
variable.  The  value  of  REFL  for  land  (0.07  ± 0.05  rms)  was 
taken  from  empirical  studies  by  Rosenkranz  (Ph.D.  Thesis, 
Department  of  Electrical  Engineering,  MIT,  1971) Hard 
limiting  was  used  to  ensure  that  no  values  of  REFL  were 
outside  the  range  0 to  1. 

The  theoretical  basis  for  calculating  the  D matrix  was 
discussed  earlier  in  this  report. 

Also,  a procedure  exists  to  estimate  the  error  which 
would  occur  for  a given  amount  of  noise  added  independently 
to  each  channel.  This  estimate  is  derived  as  follows: 

p*  = D (d  + n)  , 

where  p*  = estimated  parameter  vector, 

D = D matrix, 

d = noiseless  observation  vector, 
and  n = noise  vector. 

In  this  case,  the  expected  contribution  to  the  covariance  p* 
added  by  a non-zero  n which  is  independent  of  d is: 


Table  4.3.  Subroutines 

MH20L  - Refractive  index  of  water  from  Ray,  Applied  Optics, 

Vol . 11,  No.  8. 

Complex  refractive  index,  Temp  °K,  Freq  in  GHz, 

Salinity  fraction,  Pure  H20/Seawater/NaCl  sol.  j 

MICE  - Refractive  index  of  ice  from  Ray,  Applied  Optics, 

Vol.  11,  No.  8. 

Complex  refractive  index,  Temp  °K  (173-273), 

Freq  in  GHz  . 

SEAT  - Finds  temp,  of  sea  as  a function  of  latitude 
assuming  mid-Pacific  temperature  profile. 

Latitude,  GMT  day.  Calculated  Sea  Temperature. 

DIELEC  - Computes  complex  dielectric  constant  for  H20 
from  Stogryn's  equations  from  IEEE  Trans., 

MTT-19 , 733,  1971. 

Freq  GHz,  Temp  'K,  Salinity  fraction  by  weight 
(0.035  for  sea  water),  NaCl  sol 1 n/Seawater . 

ABH202  - Computes  H20  absorption  due  to  22.235  GHz  and 

183.3  GHz  lines;  derived  from  AB4H20  , 

R.  Pettyjohn,  Dec.  '75. 

Temp  °K , Pressure  mb,  H20  vapor  density,  Freq.  GHz. 

02ABSB  - Dry  air  absorption  coefficient  due  to  oxygen  in 
Nepers/Km  by  P.  W.  Rosenkranz,  IEEE  Trans., 

AP-24  Preprint  . 

Temp  °K,  pressure  mb,  Freq  GHz. 

RHOSAT  - Computes  saturated  Ho  vapor  density. 

Temp  ° K . 2 i 

GAMLW  - Computes  absorption  coefficient  in  Nepers/Km  due 
to  liquid  H20,  excluding  scattering. 

P.  W.  Rosenkranz  Ph.D.  Thesis  1971. 

RAIN  - Puts  cloud  models  of  Gaut  and  Reifenstein  into 
atmospheres.  Also  includes  some  interpolated 
models.  From  ERT  Technical  Report  No.  12,  Microwave 
Properties  of  Clouds  in  the  Spectral  Range  30-40  GHz, 

18  Jan  1977,  E.  C.  Reifenstein  and  N.  E.  Gaut.  ' j 
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Ef  (Dn)  (Dn)1]  = E[ (Dn)  (n^Dt) ] = E(DnntDt). 

If  the  noise  in  each  channel  is  zero  mean,  of  equal  variance, 
and  mutually  independent  then 

E (DnntDt)  = E (DaIDt)  = aDDfc, 

where  a is  a scalar  denoting  the  magnitude  of  the  noise. 

The  program  assumes  this  to  be  1,  i.e.,  the  variance  of  the 

2 

noise  is  assumed  to  be  (1°)  . Although  the  non-diagonal 

elements  are  non-zero  in  DDfc,  the  diagonal  elements  are  the 
most  interesting,  since  they  contribute  to  the  rms  error  in 
each  individual  estimate. 

4 . 6 Results  of  Computations 

The  following  parameters  were  singled  out  for  inclusion 
in  the  summary  results  presented  in  Appendix  A of  this 
report : 

Temperatures  at  1000  and  500  mb  (°K) 

Integrated  Water  Vapor  (cm) 

Integrated  Liquid  Water  (cm  x 100) 

Cloudtop  Height  (km  x 10) 

Cloudtop  Temperature  (°K) 

The  temperature  at  1000  mb  was  chosen  since  it  produces  the 
least  accurate  retrievals  of  any  pressure  level,  and  hence 
sets  a ceiling  on  performance.  500  mb  was  chosen  because 
it  has  the  most  accurate  retrievals  and  thus  establishes 
the  lowest  upper  bound  on  performance. 

Within  each  box  in  the  charts  presented  in  Appendix  A, 
the  entries  are  arranged  as  follows: 
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where,  if  is  the  parameter  to  be  estimated,  and  p^  is  the 

estimate,  p7  is  the  expected  value  of  p^,  and  E(x)  denotes 
the  expected  value  of  p. , 


A.  = (E[(p.  - Pi)2] 


',  the  rms  error, 


and 


2 ) *5 

= E [ (p^  - p^)  ] , the  a priori  variation, 

t 

= (DD  ) , the  "D  matrix  covariance," 

A 

= [ E (pi ) ] - Pi,  the  mean  error, 

E^  = p7,  the  parameter  mean. 


The  different  types  of  atmospheres  used  for  the  compari- 
sons are  best  characterized  by  the  mean  values  of  their 
parameters,  as  indicated  by  E^ . The  D matrix  is  always 
derived  from  the  #1  data  set  and  thus  the  mean  retrieval 
errors  for  the  #1  data  sets  are  generally  zero.  When 
retrievals  are  included  for  a "#2"  data  set  of  the  same 
type  as  a "#1"  data  set,  the  D matrix  calculated  from  the 
"#1"  data  set  is  applied  to  the  "#2"  data  set.  In  this 
case,  mean  errors  can  be  expected  and  the  estimates  will 
generally  not  be  so  good  as  for  the  "#1"  data  sets,  although 
the  results  can  be  taken  to  be  more  realistic. 
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4 . 7 Non-Linear  Estimation 

An  unsuccessful  attempt  was  made  at  solving  a simple 

non-linear  problem,  not  necessarily  because  it  would  be 

directly  applicable  to  this  retrieval  problem,  but  rather 

to  illustrate  the  difficulty  of  such  an  analytically-based 

retrieval  procedure.  The  problem,  simply  stated,  was  as 

follows:  A given  parameter  x has  a Gaussian  distribution 

2 2 

with  mean  y and  variance  a . (The  value  of  y and  a may  not 

necessarily  be  known.)  This  parameter  is  not  observed 

2 

directly,  but  rather  x is  observed,  corrupted  with  zero 

2 

mean  Gaussian  noise  of  variance  a . The  problem  is  to  find 
the  probability  distribution  of  this  observation.  Thus: 

2 

x is  N (y , a ) , 

2 

n is  N (0,  o ) , 

2 . 

y = x + n. 

In  general,  if  we  can  find  the  probability  distribution 

function  for  an  observation  y,  then  we  can  find,  for  exam- 

2 

pie,  the  maximum  likelihood  estimate  of  y and  a given  a 
number  of  observations;  we  could  also  find  the  least  square 
error  estimate  of  x for  each  value  of  y we  observe.  The 
problem  is  instructive  when  tackled  analytically,  since  it 
illustrates  the  intractability  of  the  non-linear  problem, 
and  since  the  simple  nature  of  this  particular  problem 
lends  itself  to  checks  via  other  methods. 

The  first  approaches  to  the  problem  involved  the  MIT 
Mathlab  computer  program  MACSYMA^Vhich  is  a program  which 
manipulates  algebraic  expressions.  It  can  perform  a large 
class  of  operations,  including  symbolic  integration  and 
differentiation,  Fourier  and  Laplace  transforms,  matrix 
operations,  etc. 

Since  y is  the  sum  of  two  random  variables  of  known 
distribution,  and  which  are  assumed  to  be  independent,  its 


\ 


M 
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probability  distribution  is  the  convolution  of  the  probabil- 
ity distributions  of  the  two  random  variables.  This  is  a 

difficult  integration  to  perform  since  x is  not  directly 
. 2 

involved  in  the  problem;  rather  x is,  and  x had  non-zero 

mean.  Originally,  the  approach  tried  was  to  take  the 

2 

Laplace  transforms  of  the  distributions  of  x and  n,  multi- 
ply them,  and  take  the  inverse  transform.  Unfortunately, 
the  inverse  transform  of  the  product  did  not  yield  to  any 
solution,  even  with  the  aid  of  MACSYMA , although  the  moment 
generating  properties  of  the  transform  did  give  the  expected 
results. 

A direct  attempt  was  then  made  to  perform  the  convolu- 
tion. Unfortunately,  due  to  a bug  in  MACSYMA  (it  cannot 
correctly  integrate  non-zero  mean  Gaussians) , the  attempt 
failed.  However,  in  about  one  day,  with  the  aid  of 
Gradshteyn  and  Ryzhik,  Tables  of  Integrals,  Series,  and 
Products , the  integration  was  performed  without  the  aid  of 
the  computer.  Unfortunately,  the  result  must  be  expressed 
in  terms  of  special  functions  (which  are  tabulated) ; it  is 
not  clear  how  this  result  could  be  used  in  estimation 
problems  as  an  analytical  function. 

Because  even  this  simple  problem  gave  rise  to  a very 
difficult  solution  of  limited  analytical  use,  it  would 
appear  that,  even  in  relatively  simple  non-linear  cases,  a 
numerical  method  of  estimation  would  seem  to  be  the  most 
productive  approach. 

4 . 8 Conclusions 

4.8.1  Temperature 

The  60  and  20  GHz  channel  set  gave  comparable  1000  mb 
temperature  retrieval  results  over  land  and  ocean  for  clear 
air  and  light  clouds,  but  was  degraded  for  heavy  clouds 
filling  the  antenna  beam.  "Heavy"  clouds  are  defined  as 
those  having,  for  a depth  of  at  least  1000  meters,  a density 
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of  at  least  1 gm/m^  of  liquid  water.  (No  clouds  with 
greater  densities  over  lesser  depths  were  included  in  the 
cloud  models  used  here.)  The  118,  100,  and  135  GHz  channel 
set  (possibly  excluding  135  GHz,  since  it  is  the  least  ef- 
fective channel  for  these  retrievals)  is  somewhat  degraded 
over  land,  but  is  slightly  improved  over  clear-air  ocean 
compared  to  the  60  and  20  GHz  set;  it  is  comparable  in  the 
light  cloud  ocean  case.  Adding  all  sets  together  (60,  20, 
118,  100,  135,  and  183  GHz  sets)  generally  offers  slight 
improvement. 

The  temperature  at  500  mb  was  tabulated  only  for  land 
retrievals.  The  tabulations  indicate  that  the  retrievals, 
for  all  channel  sets,  are  affected  by  clouds  to  a much 
lesser  degree  than  the  1000  mb  temperatures,  which  is  as 
expected.  It  also  indicates  that  the  channel  sets  which 
include  both  60  GHz  and  channels  above  60  GHz  are  somewhat 
more  immune  to  errors  caused  by  clouds  than  is  the  60  GHz 
set  alone. 

4.8.2  Water  Vapor 

As  expected,  water  vapor  retrievals  were  poor  over  land. 
Over  ocean,  excellent  results  were  obtained  with  the  60  and 
20  GHz  channel  set  over  clear  air,  with  performance  de- 
graded somewhat  by  light  clouds  and  very  considerably  by 
heavier  clouds.  The  118,  100,  and  135  GHz  channel  set 
yielded  performance  over  clear-air  ocean  only  slightly 
worse  than  the  60  and  20  GHz  set;  the  183  GHz  set  did  not 
come  into  its  own  until  clouds  were  added;  the  118,  100, 

183,  and  135  GHz  set  offered  comparable  performance  to  the 
60  and  20  GHz  set  except  that  it  was  more  severely  de- 
graded by  light  clouds.  The  complete  28  channel  set  pro- 
vided little  improvement  over  the  60  and  20  GHz  set  alone. 


4.8.3  Liquid  Water 

Liquid  water  retrievals  seem  to  be  possible  with  this 
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simple  linear  retrieval  scheme  only  over  ocean  and  only  for 
light  clouds.  The  results  suggest  that,  in  this  case,  the 
60  and  20  GHz  channel  set  performs  slightly  better  than  the 
183,  118,  100,  and  135  GHz  channel  set.  The  results  also 
suggest  that  there  are  strong  non-linearities  which  make 
the  simple  retrieval  system  used  here  less  appropriate  than 
it  was  for  the  other  parameters.  Future  development  of  non- 
linear retrieval  schemes  may  yield  significant  improvement 
in  these  estimates  of  liquid  water. 

4.8.4  Cloudtop  Height  and  Temperature  at  the  Cloudtop 

In  this  analysis  the  cloudtop  height  was  defined  to  be 
either  the  ground,  or  the  highest  altitude  at  which  at  least 
0.3  g/m"^  of  liquid  water  exists,  whichever  is  higher.  Over 
land  the  cloudtop  height  estimates  are  not  very  impressive 
in  the  light  cloud  case;  for  the  heavier  cloud  cases,  it  is 
found  that  in  the  tropics  the  higher  frequencies  do  not  do 
so  well  as  the  lower  frequencies,  whereas  at  mid-latitude, 
the  reverse  is  true.  Over  ocean,  the  retrievals  are  much 
better,  especially  for  the  light  cloud  cases.  For  tempera- 
ture at  cloudtop,  the  results  follow  roughly  the  same  pat- 
tern. For  heavy  clouds,  the  retrievals  are  remarkably  good. 
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Appendix  A 

Tabulated  Results:  Retrieval  Accuracies 


Midlatitudes : 

Tropics : 

1000  mb  Temperature  (Over  Ocean) 

1000  mb  Temperature  (Over  Ocean) 

Midlatitudes : 

Tropics : 

Temperature  at  Cloudtop  (Over  Ocean) 

Temperature  at  Cloudtop  (Over  Ocean) 

Midlatitudes : 

Tropics : 

Water  Vapor  (Over  Ocean) 

Water  Vapor  (Over  Ocean) 

Midlatitudes : 

Tropics : 

Liquid  Water  (Over  Ocean) 

Liquid  Water  (Over  Ocean) 

Midlatitudes : 

Tropics : 

Cloudtop  Height  (km  x 10)  (Over  Ocean) 

Cloudtop  Height  (km  x 10)  (Over  Ocean) 

Midlatitudes : 

Tropics : 

500  mb  Temperature  (Over  Land) 

500  mb  Temperature  (Over  Land) 

Midlatitudes : 

Tropics : 

1000  mb  Temperature  (Over  Land) 

1000  mb  Temperature  (Over  Land) 

Midlatitudes : 

Tropics : 

Integrated  Water  Vapor  (Over  Land) 

Integrated  Water  Vanor  (Over  Land) 

Midlatitudes : 

Tropics : 

Cloudtop  Height  (Over  Land) 

Cloudtop  Height  (Over  Land) 

Midlatitudes : 

Tropics : 

Cloudtop  Temperature  (Over  Land) 

Cloudtop  Temperature  (Over  Land) 
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Midlatitudes:  1000  mb  Temperature  (Over  Ocean) 
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Tropics:  1000  mb  Temperature  (Over  Ocean) 


N 

X 

r- 


oo  in  o 

O 

<N 

co  in  cn 

*"H  | 

CO 

O 

O 

rH  in  ID 

o 

CN 

CO  ID 

CO 

rH  CO 

H \ • 

• 

• 

\ • • 

\ 

• 

• 

• \ • 

• 

• 1 

S • 

• 

• \ 

rH  in  rH 

o 

VD 

O rH  o 

00 

rH 

rH 

o 

00  rH 

o 

ID 

CN  rH 

o 

00  CO 

• 

on 

• 

ON 

• 

ON  • 

on 

• 

on  • 

VI  <N 

CN 

cn 

CN 

CN 

CN  CO 

CN 

CO 

CN  CN 

CO 

00  rH 

CN  (1) 

rH 

rH 

rH 

r- 

rH 

C 

• 

• 

• 

• 

• 

rH  C 

in  in 

O 

CN 

co  in 

rH 

rH 

CO  rH 

O 

rH  in  CN 

o 

CN 

CO  CN 

rH 

rH  rfl 

\ • 

• 

• 

\ • 

• 

\ • 

• 

• \ • 

• 

• i 

\ . 

• 

• 

< -C 

CN  rH 

O 

ID 

O'*  rH 

O 

00 

O rH 

o 

oo  in  (N 

o 

ON  CN 

O 

00 

u 

• 

ON 

• 

on  ; 

• 

ON  • 

ON 

• 

ON 

CN 

CN 

rH 

CN 

CN 

CN  CO 

CN 

CN 

CN 

r"* 

in  oo  o cn 


rH 

rH 

rH 

• 

• 

• 

• 

CN  CO  00 

00 

rH 

co 

rH 

O 

rH  in 

rr 

o 

CN 

CO  Tf 

CO 

rH 

• \ • 

• 

• 

\ 

• 

• 

• \ 

• 

• 

\ . 

• 

• 

ID  ON  rH 

o 

00 

CN 

rH 

o 

00  00 

rH 

o 

VD  ! 

CN  rH 

o 

00 

ON 

ON 

■ 

ON 

ON  ; 

• 

ON 

CN  CN 

CN 

CN 

CN  CO 

CN 

CO 

CN 

rH 

rH 

r- 

rH 

• 

• 

• 

• 

CN  CO  CN 

CN 

rH 

CO  CO 

O 

rH  in 

O 

CN 

co  r- 

rH 

• \ • 

• 

• 

\ • 

• 

• \ 

• 

• 

• 

\ • 

• 

• 

VD  ON  rH 

O 

00 

CO  rH 

o 

00  rH 

CN 

o 

VD 

rH  CN 

o 

00 

ON  • 

ON 

• 

ON 

ON 

• 

ON 

CN  rH 

CN 

CN 

CN  •O' 

CN 

co 

CN 

4->  CO 

in 

rH 

O 

CN 

CO 

rH 

CN 

rH 

co 

CN 

o 

rH  in 

O 

o 

CN 

co 

o 

CO 

rH 

* *H  rH 

\ 

• 

• 

• 

\ 

• 

• 

• 

• 

• 

• \ 

• 

• 

• 

\ 

• 

• 

• 

CO  £ 

co 

rH 

o 

VD 

00 

rH 

O 

00 

CO 

rH 

o 

00  rH 

fSJ 

o 

VD 

CN 

CN 

o 

00 

oo  O i 

• 

ON 

• 

ON 

• 

ON  • 

ON 

• 

ON 

rH 

. CN 

CN 

1 r“^ 

CN 

1 

CN  -O' 

CN 

CO 

CN 

rH 

m 

CN 

O 

CN 

CO 

CN 

CN 

rH 

• 

\ 

• 

• 

• 

\ 

• 

• 

• 

00 

o 

CN 

o 

VD 

rH 

CN 

O 

00 

ON 

• 

ON 

• 

ON 

CN 

CN 

CO 

CN 

rH  in  o 

O 

CN 

CO 

o 

CO 

rH 

• \ • 

• 

• 

\ 

• 

• 

• 

00  rH  CN 

o 

VD 

rH 

CN 

o 

CO 

ON 

ON 

• 

ON 

CN  ^ 

CN 

CO 

CN 

in 

O 

<N 

co 

rr 

CN 

rH 

CO 

co 

o 

rH  m CO 

o 

CN 

CO  co 

CN 

rH 

«*  * 

\ • 

• 

• 

• 

• 

• 

• 

• 

• \ • 

• 

• 

1 ^ * 

• 

• 

CO  O 

CO  rH 

o 

VD 

ON 

rH 

O 

00 

CN 

rH 

o 

00  O CN 

o 

VD 

O CN 

o 

00 

00  o 

• 

ON 

• 

ON 

• 

ON 

ON 

• 

ON 

rH  rH 

CN 

CN 

rH 

CN 

CN 

cn 

CN 

CO 

CN 

co  o 
\ • • 
rH  rH  O 

(N 


>,  c 

>1  c 

to 

> -H 

> *H 

TJ 

CO  «0 

<0  <0 

w 

3 

0)  X 

ai  « 

-U 

(0 

in 

0 

X 

X 

W W 

TJ 

Tl 

c« 

m 

u u 

3 

3 

u 

LB 

0 0 

0 

0 

If)  rH 

in  in 

A 0 

rH  rH 

<-!  IN 

•U 

^ TJ  ♦ 

•M  T3  ♦ 

T3  a 

u ♦ 

U =* 

■C 

x:  a 

X 3 

a c 

O’ 

On  0 

O'  0 

0 * 

0 

0 

•H 

•r|  rH 

•H  rH 

rH  0 

z 

z 

>3 

*1  u 

J u 

U Q 

298.1  298 


Tropics:  Temperature  at  Cloudtop  (Over  Ocean) 
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